General Dynamics of Topology and Traffic on Weighted Technological Networks 
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For most technical networks, the interplay of dynamics, traffic and topology is assumed crucial to 
their evolution. In this paper, we propose a traffic-driven evolution model of weighted technological 
networks. By introducing a general strength-coupling mechanism under which the traffic and topol- 
ogy mutually interact, the model gives power-law distributions of degree, weight and strength, as 
confirmed in many real networks. Particularly, depending on a parameter W that controls the total 
weight growth of the system, the nontrivial clustering coefficient C, degree assortativity coefficient 
r and degree-strength correlation are all in consistence with empirical evidences. 
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The past few years have witnessed a great devotion of 
physicists to understand and characterize the underly- 
ing mechanisms of complex networks including the In- 
ternet jj, the WWW 0, the scientific collaboration 
networks fSCN) H and world-wide airport networks 
So far, researches on networks have mainly 
focused on unweighted graphs. Barabasi and Albert have 
proposed a well-known model (BA model) that intro- 
duces the degree preferential attachment mechanism to 
mimic unweighted growing networks |3]. Most recently, 
the availability of more complete empirical data has al- 
lowed scientists to consider the variation of the weights 
of links that reflect the physical characteristics of many 
real networks. Obviously, there is a need for a mod- 
elling approach to complex networks that goes beyond 
the purely topological point of view. Alain Barrat, et 
al. presented a model (BBV model) that integrates the 
topology and weight dynamical evolution to study the 
growth of weighted networks Their model yields 

scale-free properties of the degree, weight and strength 
distributions, controlled by an introduced parameter 5. 
However, its weight dynamical evolution is triggered only 
by newly added vertices, resulting in few satisfying in- 
terpretions to the collaboration networks or the airport 
systems. In fact, the dynamics and properties of social 
and technological networks are quite different and should 
be addressed individually. It is well-known that networks 
are not only specified by their topology but also by the 
dynamics of weight (e.g. information flow) taking place 
along the links. For instance, the heterogeneity in the 
intensity of connections may be very important in un- 
derstanding technological systems. The amount of traf- 
fic characterizing the connections of communication sys- 
tems or large transport infrastructure is fundamental for 
a full description of these networks Take the WAN 
for example: each given edge weight Wij (traffic) is the 
number of available seats on direct flight connections be- 



tween the airports i and j. Weighted networks are often 
described by an adjacency matrix Wij which represents 
the weight on the edge connecting vertices i and j, with 
i,i = 1, . . . ,N, where N is the size of the network. We 
will only consider undirected graphs, where the weights 



are symmetric {wi 



As confirmed by measure- 



ments, complex networks often exhibit a scale- free de- 
gree distribution P{k) ~ fc-^ with 2< 7 <3 The 
weight distribution P{w) that any given edge has weight 
w is another significant characterization of weighted net- 
works, and it is found to be heavy tailed, spanning sev- 
eral orders of magnitude |lO| |. A natural generalization 
of connectivity in the case of weighted networks is the 
vertex strength described as = ^ 



Wi 



where the 



sum runs over the set r(i) of neighbors of node i. The 
strength of a vertex integrates the information about its 
connectivity and the weights of its links. For instance, 
the strength in WAN provides the actual traffic going 
through a vertex and is obvious measure of the size and 
importance of each airport. Empirical evidence indicates 
that in most cases the strength distribution has a fat 
tail @, similar to the power law of degree distribution. 
Highly correlated with the degree, the strength usually 
displays scale- free property s ^ pdl IT^ . 

The previous models of complex networks always in- 
corporate the (degree or strength) preferential attach- 
ment mechanism, which may result in scale-free proper- 
ties. Essentially, this mechanism just describes interac- 
tions between the newly-added node and the old ones. 
Actually, such interactions also exist between old nodes. 
Perhaps, the most reasonable and simplest way to ex- 
press such interactions is by the product form of related 
vertex strengths, i.e. the pairwise interaction between 
vertices i and j is proportional to SiSj (strength coupling 
form). Let's review BA model: a new vertex n is added 
with m edges that are randomly attached to an existing 
vertex i according to the degree preferential probability, 
which can be written in the product form of degrees 
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Analogously, in BBV networks one can rewrite the 
strength preferential probability: 



n„ 



(2) 



We argue that such interactions (actually driven by traf- 
fic) exist between old vertices in the same way, and will 
considerably affect the flows between them: First, new 
edges should be allowed to add between old nodes; sec- 
ond, the pre-existing traffic flows along the links will be 
updated with the growth of networks. Indeed, physi- 
cal interaction of nodes plays a crucial role in determin- 
ing the network topology during its dynamical evolution. 
Our above perspectives have been partly inspired by the 
work of Dorogovtsev and Mendes (DM) who pro- 
posed a class of undirected and unweighted models where 
new edges are added between old sites (internal edges) 
and existing edges can be removed (edge removal). 

In the letter, we present a model for weighted techno- 
logical networks that considers the topological evolution 
under the general trafflc-driven interactions of vertices. It 
can mimic the reinforcement of internal connections and 
the evolution of many infrastructure networks. The di- 
versity of scale-free characteristics, nontrivial clustering 
coefficient, assortativity coefficient and strength-degree 
correlation that have been empirically observed can be 
well explained by our microscopic mechanisms. More- 
over, in contrast with previous models where weights are 
assigned statically [l4L Il5l | or rearranged locally (BBV 
model), we allow the flows to be widely updated. 

The model starts from an initial configuration of A^o 
vertices connected by links with assigned weight wq ■ The 
model is defined on two coupled mechanisms: the topo- 
logical growth and the strengths' dynamics: 

(i) Strengths ' Dynamics. From the beginning of the 
evolution, all the possible (existing or not) connections 
at each time step are supposed to update their weights 
according to the strength-coupling mechanism: 



"I] 



1, with probability Wpij 
, with probability 1 — Wpij , 



where 



(3) 



(4) 



integrates the strength coupling of vertices i and j, and 
determines the increment probability of weight Wij (if i 
and j are unconnected, Wij — 0). The total weight of 
the edges in statistical sense is modified by the amount 
< X]i<j ^'^ij which is assumed constant for sim- 

plicity. This parameter reflects the growing speed of the 
network's total traffic load; for instance, the increasing 
rate of total information flow in a communication sys- 
tem. The always growing traffic plays the driving role 
in network evolution. One may notice that Wpij is very 
likely to exceed one if the initial number of nodes A^o is 
small. When Wpij exceeds one, it is automatically as- 
sumed to be one. This treatment of Wpij will probably 




FIG. 1: Illustration of the evolution dynamics. A new 
node n connects to a node i with probability proportional 



to Si/E, 



The thickness of nodes and links respectively 



represents the magnitude of strength and weight. New con- 
nections (dashed) can be built between pre-existing nodes and 
the bilateral links represent the traffic growing process along 
links under the general mechanism of strength couplings. 



affect the initial network evolution, while it is not signif- 
icant for discussing the statistical measures, as they are 
almost independent of initial states. 

(ii) Topological Growth. At the same time step, a new 
vertex n is then added with m edges that are randomly 
attached to an existing vertex i according to the strength 
preferential probability Iln^i (Eq- (2)). The weight of 
each new edge is also fixed to wq. In fact, the strength 
preferential attachemnt is essentially the same with the 
mechanism traffic-driven interactions we have argued. 

The network provides the substrate on which numer- 
ous dynamical processes occur. In previous models, traf- 
fic was often assumed just as an appendix to the network 
structure. Actually, traffic and the underlying topology 
are mutually correlated and it is very important to define 
appropriate quantities and measures capable of captur- 
ing how all these ingredients participate in the formation 
of complex networks 9]. Technology networks provide 
large empirical database that simultaneously captures 
the topology and the dynamics taking place on it. For 
Internet, the information traffic between routers (nodes) 
can be represented by the corresponding edge weight. 
The total traffic that each router deals with can be de- 
noted by the node strength, which also represents the 
importance of given router. The increasing information 
fiow as an internal demand always spurs the expansion 
of technological networks. Specifically, the largest contri- 
bution to the growth is given by the emergence of links 
between already existing nodes. This clearly points out 
that the Internet growth is strongly driven by the need 
of a redundancy wiring and an increasing need of avail- 
able bandwidth for data transmission 12]. On one end, 
newly-built links (between existing routers) are supposed 
to preferentially connect high strength routers, because 
otherwise, it would lead to the unnecessary traffic conges- 
tion along indirect paths that connect those high strength 
nodes. Naturally, traffic along existing links between high 
strength routers, in general, increases faster than that be- 
tween low strength routers. All the points are reflected 
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in our strength-coupling mechanism. On the other end, 
new routers preferentially connect to routers with larger 
bandwidth and traffic handling capabilities (the strength 
driven attachment). Those phenomena also exist in air- 
port system, power grid, and railroad network; and they 
could be explained by the traffic-driven mechanism of our 
model. For power grid and railroad network, the cost by 
distance has a distinct effect to their topological prop- 
erties. Their degree distributions, for example, are not 
scale-free. In a word, topology and traffic interact with 
each other in networks under general interactions of ver- 
tices driven by the internal increasing traffic demand. 

The model time is measured with respect to the num- 
ber of nodes added to the graph, i.e. t ^ N — Nq, and 
the natural time scale of the model dynamics is the net- 
work size N. In response to the demand of increasing 
traffic, the system must expand. With a certain size, one 
technological network assumably has a certain ability to 
handle certain traffic load. Therefore, it could be rea- 
sonable to suppose that the total weight on the networks 
increases synchronously by the natural time scale. That 
is why we assume as a constant. This assumption also 
bring us the convenience of analytical discussion. Using 
the continuous approximation, we can treat fc, w, s and 
the time t as continuous variables [10. Then Eq. (3) 
indicates: 

dwij 2WsiSj 2WsiSj 

The strength s; of vertex i can increase either if a new- 
added node connects to i by the topological growth dy- 
namics or any possible (exsiting or not) connections to i 
are updated by the strengths' dynamics: 

dsj _ J2]{^i) '^WsjSj rnsj _ 2W + m Sj 

where the last expressions recovered by noticing that 
J2iSi{t) ~ 2{m + W)t. From Eq. (6), one can ana- 
lytically obtain the power-law distribution of strength 
(P(s) s~") with the exponent |l|i|: a = 2 + m/{m + 
2W). Obviously, when W = Q the model is topologi- 
cally equivalent to the BA network and the value a — i 
is recovered. For larger values of VF, the distribution is 
gradually broader with a — > 2 when W oo. 

We performed numerical simulations of networks gen- 
erated by choosing different values of W and fixing Nq — 
3, m = 3 and — 1. We have checked that the scale- 
free properties of our model networks are almost inde- 
pendent of the initial conditions. Numerical simulations 
are consistent with our theoretical predictions, verifying 
again the reliability of our present results. Fig. 2(a) gives 
the probability distribution P{s) ~ s", which is in good 
agreement with the theoretical predictions. We also re- 
port the average strength Si of vertices with degree ki, 
which displays a nontrivial power-law behavior s ~ fc'^ as 
confirmed by empirical measurement. Unlike BBV net- 
works (where /? = 1), the exponent (3 here varies with 




s k. 



FIG. 2: (a) Probability distribution P{s). Data are consis- 
tent with a power-law behavior s~" . In the inset we give the 
value of ct obtained by data fitting (filled circles), together 
with the analytical expression a = 2 -\- mjim -|- 2M^)(line). 
The data are averaged over 20 networks of size N=5000. (b) 
Strength Si versus ki for different W (log- log scale). Linear 
data fitting gives slope 1.04, 1.17, 1.25 and 1.30 (from bottom 
to top), demonstrating the correlation of s ~ fc''. 
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FIG. 3: (a) Probability distribution of the degrees P(fc) ~ 
k~^ . (b) Probability distribution of the weights P(io) ~ w~'^ . 
The data are averaged over 20 networks of size N=5000. 

the parameter TF' in a nontrivial way as shown in Fig. 
2(b). Moreover, the major difference between our model 
and DM network is reflected in the strength-degree cor- 
relation graph. Though DM model allows the emergence 
of internal edges, it could not mimic the reinforcement 
of pre-existing connections in that it is unweighted. The 
nontrivial s ^ k^ correlation demonstrates the signifi- 
cant part of weight increment along existing edges, and 
thus implies that our model is reasonable in the light 
of traffic flow. More importantly, one could check the 
scale-free property of degree distribution {P{k) ^ k~'^) 
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FIG. 4: (a) Clustering coefficient C depending on the pa- 
rameter W . In the inset, we report the evolution of clustering 
coefficient (or C versus A'^) which converges soon, (b) Degree- 
degree correlation r depending onW. In the inset, we report 
its evolution which converges soon. 



by combining s ^ with P{s) ^ s^". Considering 
P{k)dk = P{s)ds, the exponent 7 is easily calculated: 
7 = f3{a — 1) -I- 1. The scale- free properties of weight 
and degree obtained from simulations are presented in 
Fig. 3(a) (b). Finally, it is worth remarking that for the 
BA networks, the clustering coefficient is nearly zero, far 
from the practical nets that exhibit a variety of small- 
world properties. In the present model, however, clus- 
tering coefficient C is found to be a function of W (Fig. 
4(a)), also supported by empirical data of a broad range. 

In the social networks, connections between people 
may be assortative by language or by race. Mixing can 
also be disassortative, i.e. vertices in the network prefer- 
entially form connections to others unlike them. Newman 
proposed some simple measures for these types of mix- 
ing, which we call assortativity coefficients jT^- the 
case of mixing by vertex degree, a remarkable pattern 



emerges. Almost all the social networks studied show 
positive assortativity coefficients while all others, includ- 
ing technological and biological networks, show negative 
coefficients. It is not clear if this is a universal prop- 
erty; the origin of this difference is not understood ei- 
ther. In our views, it represents a feature that should be 
addressed in each network individually. We argue that 
the adaptive evolution of topology in response to the in- 
creasing traffic is the major cause of disassortative mixing 
of technological networks. Using the formula defined in 
Ref. (Eq. (26)), we calculate the degree assortativity 
coefficient (or degree-degree correlation) r of the graphs 
generated by our model. Simulations given in the Fig. 
4(b) are supported by empirical measurements The 
restriction of our model to technological networks is for 
that there are few empirical data for statistical analysis 
on "weighted" biological networks, where many interact- 
ing mechanisms are far from present knowledge either. 
But hopefully, our model may be very beneficial for fu- 
ture understanding or characterizing biological networks 
and social ones, as it generates many topological prop- 
erties observed in those real networks. Due to its ap- 
parent simplicity and the variety of controllable results, 
we believe that some of its extensions will probably help 
address the other two classes of networks. 

In conclusion, the universal interactions of nodes and 
internal traffic demand of the system will determine the 
topology evolution of technological network. This gen- 
eral traffic-driven mechanism provides a wide variety of 
scale-free behaviors, clustering coefficient and nontrivial 
correlations, depending on the parameter W that governs 
the total weight growth. All the results are supported by 
empirical data. Therefore, our present model for all prac- 
tical purposes will demonstrate its applications in future 
weighted network research. 
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